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Abstract 

Computational fluid dynamics is used to simulate the wake flow and afterbody 
heating of the Project Fire II ballistic reentry to Earth at 11.4 km/ sec. Laminar 
results are obtained over a portion of the trajectory between the initial heat pulse 
and peak afterbody heating. Although forebody convective heating results are in 
excellent agreement with previous computations, initial non-catalytic predictions of 
afterbody heating were about a factor of two below the experimental values. How- 
ever, analogy with currently produced thermal protection materials suggests that 
significant catalysis may be occurring on the afterbody heat shield. Computations 
including finite-rate catalysis on the afterbody surface are in good agreement with 
the data over the early portion of the trajectory, but are conservative near the 
peak afterbody heating point, especially on the rear portion of the conical frustum. 
Further analysis of the flight data from Fire II shows that peak afterbody heat- 
ing on the frustum occurs before peak forebody heating, a result that contradicts 
computations and flight data from other entry vehicles. This result suggests that 
another mechanism, possibly pyrolysis, may be occurring during the later portion 
of the trajectory, resulting in less total heat transfer than the current predictions. 


Int reduction 

Thirty six years later, the Project Fire 11 ballistic 
reentry to Earth at a nominal velocity of 11 A krn/s 
remains one of the best sources of aerothermal heat- 
ing data for the design of sample return capsules. The 
data from this flight experiment 1 "’ 3 encompass both 
the I hermochemical non equilibrium and equilibrium 
flow regimes and include measurements of both radia- 
tive and total heating on the forebody and afterbody. 
Because of the quality of this data, a number of re- 
searchers have performed computational fluid dynam- 
ics (CFO) simulations of the forebody of the Fire II 
entry vehicle, 4-8 with generally good results. In par- 
ticular, Olynick et.al. 4 coupled a Navier-Stokes solver 
with a nonequilibrium radiation code and showed good 
agreement for stagnation point heat transfer over the 
Fire 11 Irajectory between 1 God and 1651 seconds (77 
km to :i7 km), a period that encompassed the peak 
fore body heating point (/. = 1645*). However, in most 
cases the primary motivation of the previous work 
was to model the coupling between shock layer ra- 
diation and aerotherrnodynamics, and thus the sim- 
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illations concentrated on the forebody flow only. To 
our knowledge there have been no prior published at- 
tempts to reproduce the afterbody heating data pre- 
sented in Ref. 5. 

Uncertainties in afterbody heating predictions can 
have a major impact on Thermal Protection System 
(TPS) material selection and weight. Conservatism in 
the afterbody heat shield design will shift 1 ho center 
of gravity backwards, reducing stability and in some 
cases necessitating ballast in the nose. Current design 
practices for afterbody heatshields typically assume 
a laminar, fully-catalytic, non ablating surface. The 
predictions thus obtained are augmented by a large 
factor of safety to account for turbulent transition and 
uncertainty in the baseline computations. The main 
reason for this uncertainty is a sparsity of data for val 
idation of our computational tools, (f round test data 
are typically complicated by sting interference effects, 
although a limited amount of shock tunnel data are 
available.' Little flight data at the appropriate entry 
velocities exist, and recent attempts to propose dedi 
cated flight experiments have failed to reach fruition. 
Therefore, it is important to understand the limited 
flight data that are available to improve the design fi 
delity of the next generation of Earth and planetary 
entry vehicles and to assess the need for additional 
flight data. In this paper we take a first step towards 
an understanding of the afterbody heating data from 
the Project Fire II flight experiment. 
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Figure 1. Schematic of the Fire II reentry vehicle showing multiple forebody heatshields and afterbody instru 
merit placement. Station numbers correspond to those in Table 1. (Adapted from Ref. 3) 


Geometry and Instr umen tati on 

The Fire II reentry vehicle consisted of a multi- 
layer configuration made up of three phenolic asbestos 
heat shields sandwiched between beryllium calorime- 
ters. Figure 1, taken from Ref. 3, is a schematic of the 
vehicle showing the configuration of each calorimeter. 
The first two calorimeters and their associated heat 
shields were designed to be ejected after the onset of 
melting, yielding three separate data-gathering peri 
ods. Forebody instrumentation consisted of thermo 
couples and forward looking radiometers. (See Refs. 1- 
2 for details on forebody instrumentation and place- 
ment.) 

The (16° included angle conical afterbody section 
was constructed of a fiberglass shell supporting a layer 
of phenolic-asbestos heat protection material. A thin 
surface coating of silicon elastomer was added for pre- 
launch moisture protection. The conical frustum por- 
tion of the afterbody was instrumented with a sym- 
metrical array of 12 gold calorimeters, distributed at 
three circumferential locations and four ^-stations on 
the frustum, as shown in Fig. 1. In addition, a single 
rearward-facing radiometer was placed on the frustum, 
also shown in Fig. 1. Details of the afterbody instru 
inenl placement are given in Table 1. The published 
uncertainty of the thermocouples at the base of each 
calorimeter was ±28 K over the measured tempera- 
ture range. 3 Raw temperature data were converted to 
estimates of total surface heat transfer by a curve fit- 


ting procedure. 10 Uncertainty in the final heat transfer 
obtained is not given, however we assume here that 
the uncertainty in the computation of the tempera- 
ture gradient is of the same order as the uncertainty 
in T. This, coupled with published uncertainties in 
other steps of the data reduction process, 10 yield an 
estimated total uncertainty in heat transfer of ±15%. 
The rearward facing radiometer indicated zero radia- 
tion during the entire entry with the exception of one 
brief pulse correlated to a forebody heatshield ejection 
event. Since the intensity threshold of the radiometer 
was approximately 1 W/cirr in the 0.2 to TO micron 
range, and since the single flash indicated that the ra- 
diometer was indeed functioning, it is concluded that 
the afterbody heating rates consisted mainly of convec- 
tive heating. 3 Further details on the instrumentation 
and data reduction procedures are available in Refs. 3 
and 10. 

Although the published reports give sufficient infor- 
mation to accurately determine the forebody geome- 
try, complete information for the entire vehicle is more 
difficult to obtain. By compiling geometrical data from 
several sources 1,3,11 we were able to construct an ac- 
curate (although somewhat simplified) outer mold line 
(OML) for the first heat shield configuration, shown 
in Fig. 2. The primary simplification that was made 
during construction of the OML was the elimination 
of the C-band antenna, which can be seen at the base 
of the vehicle in Fig. 1. This modification was made to 


i 



AlAA 2001-3065 3 


Station 

Sensor Type 

x/L 

2 

Beginning of Frustum 

0 

* i 
>> 

;t Calorimeters 
4> = 0°, 120°, 240° 

0.19 

4 

1 Radiometer 
6 = 203.45° 

0.32 

5 

3 Calorimeters 
<$ = 0°, 120°, 240° 

0.38 

t; 

3 Calorimeters 
<j> = 0°, 120°, 240° 

0.56 

7 

3 Calorimeters 
= 0°, 120° ,240° 

1 Pressure Sensor, (j> = 265° 

0.70 

8 

Theoretical cone apex 

1.0 


Table 1 Instrument placement on Fire II afterbody. 
Station numbers and xjL locations refer to Fig. 1. <j> 
refers to the circumferential location of the sensor. 


allow for axisymmetric solutions to be obtained (the 
antenna is three-dimensional), and also because exact 
dimensions for the antenna could not be located. Since 
all instrumentation was placed on the conical frustum 
well away from the C band antenna, this approxima- 
tion should not greatly affect the- computed results. In 
addition, a small backward facing step resulting from 
the attachment of the third forebody heatshield to the 
vehicle was eliminated, again due to a lack of geomet- 
ric information. We are still working to completely 
define these details of the OML. 

As stated above, 1 he forebody consisted of three 
calon meter/heat shield combinations that were de- 
signed to be ejected after the onset of melting. This 
configuration effectively divided the data-gathering 
window into three discrete periods, interspersed with 
brief periods where the heatshield ejection process 
and resultant flowfield interaction made accurate data 
gathering impossible. The data periods bracket the 
peak forebody heating porlion of the entry trajectory. 
Although the nominal angle of attack of the vehicle 
was 0°, the heatshield ejection process introduced sig- 
nificant off-axis forces to the reentry vehicle, resulting 
in successively larger angles of attack with time. Based 
on flight telemetry data, the angle of attack was less 
than 1° for the first data period, 5° during the second, 
and as much as 11° fluring the third period. J For this 
and other reasons (discussed below), the analysis pre- 
sented here is limited to the first data period, which 
spans the onset of significant afterbody heating to the 
peak heating period on the conical frustum. 



Figure 2. Relevant dimensions of the Fire 11 reentry 
vehicle outer mold line (OML) as modeled in the CFD 
simulations. All dimensions are in centimeters. 

Physical and Numerical Modeling 

A detailed description of the relevant equations 
for reentry flows has been presented in a number of 
sources 12-14 and will not be repeated here. Velocities 
during the portion of the Fire II trajectory considered 
here were in excess of 1 1 km/sec at an altitude range 
between 75 and GO km. At these conditions the post- 
shock flowfield is in a state of thermal and chemical 
nonequilibrium, and a significant amount of ionization 
is expected. Therefore, the flow is modeled using an 
11-species (A^, Oy , NO, /V, O, A5>+, 0‘j+, NO+ , N+, 
0+, e) 47 reaction finite-rate air chemistry model. The 
chemical source terms are modeled using rates from 
Park, 15 with the exception of the rates for the elec- 
tron impact ionization reactions, which are taken from 
Wilson. 16 The governing temperature for the electron 
impact ionization reactions is assumed to be the heavy 
particle translational temperature (T) based on earlier 
work by Olynick. 4 The flow is modeled assuming ther 
mal non equilibrium, with separate equations for con- 
servation of vibrational and total energies. The cou- 
pling between the translational and vibrational energy 
modes is modeled using a Landau Teller formulation, 
where relaxation times are obtained from Millikan and 
White, 17 assuming simple harmonic oscillators. 

Viscous transport and thermal conductivity are 
modeled using the collision cross-section fits of Gupta. 18 
The bifurcation method 19 is used to model binary dif- 
fusion, with the reference binary diffusion coefficient 
computed using the Gupta fits. This method 
allows for variations in species diffusion coefficients to 
be modeled without sacrificing the requirement that 
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the diffusion velocities sum identically to zero. The 
diffusion velocity of the electrons is computed from 
the requirement of zero current density in the flow. 
The entire flowfield is assumed to be laminar, which 
should be a good assumption, given the low freest ream 
Reynolds numbers (see Table 2) during the early por 
tion of the entry. The assumption of laminar flow in 
the near wake was validated using blunt body separa- 
tion shear layer and inner wake transition correlations 
given by Lees. 20 Even for the final (highest Re) trajec- 
tory point simulated the local Reynolds number in the 
shear layer (/?e/) was more than an order of magnitude 
below the transition Reynolds number, and the wake 
transition correlation indicated laminar flow for more 
than eleven body diameters downstream. Since the 
sonic point in the neck region was only about 2.5 di- 
ameters downstream, we conclude that transition will 
not affect the separated flow region. 

Surface catalysis is modeled using a diffusion lim- 
ited approach, 21 in which the rate of production of 
a species at the wall is balanced by the rate of mass 
diffusion of the constituent species, ie. 


where (>, is the species density, v a \ w is body-normal 
diffusion velocity of species s at the wall and w 9 
is the rate of production, modeled with a first order 
reaction rate as 



where 7£ is the universal gas constant, T w is the wall 
temperature, M., is the species molecular weight, and 
7 ,, is an accommodation coefficient representing the 
fraction of atoms or ions that recombine upon reaching 
the wall. 

We assume that the surface of the beryllium fore- 
body is completely oxidized by the high concentration 
of O atoms in the shock layer, and is thus non-catalytic 
to neutral species ( 7 „ fll/ra ^, = 0). However, the sur- 
face is assumed to be uncharged, and thus fully cat- 
alytic to ion recombination ( 7 j on ,, — l)- 21 This as- 
sumption was shown previously 4 to yield stagnation 
point heat transfer predictions in good agreement with 
the flight data. In addition, while the catalycity of 
oxidized beryllium in a reentry environment is not 
known, there is a large amount of experimental data 
indicating that oxidized metals are much less catalytic 
than their “clean’ 1 counterparts, 22-2,1 and that an ox- 
ide with low electrical conductance, such as beryllium 
oxide, should be a poor catalyst. The afterbody sur 



b) 

Figure 3. Grid topologies used in this work, enlarged 
to show detail near the body surface, a) Single block 
hyperbolic grid and b) Five block grid with compact 
enrichment in wake region. White lines delineate ele- 
mental blocks. 

face coating is silicon elastomer, which was also ini- 
tially assumed to be non-catalytic to neutral species 
and fully catalytic to ions. This assumption will be 
examined further in the results section. An isother- 
mal wall boundary condition was applied on the highly 
conductive metallic forebody, with the wall tempera- 
ture at each time period taken from the thermocouple 
data. Both a cold wall and a radiative equilibrium 
wall boundary condition were considered for the after 
body, with details given in the results. A no-slip wall 
velocity boundary condition was applied on the entire 
surface. Freestream conditions were based on atmo- 
spheric soundings taken immediately after the flight. 11 

The CFD code used to obtain solutions for this pa- 
per was DPLR-Axi. 24 DPLR-Axi is a parallel multi- 
block axisymmetric extended Navier Stokes .solver, 
based on algorithms developed in GIANTS, 4 that in- 
corporates all of the physical models discussed above. 
The Euler fluxes are computed using a modified form 
of Steger-Warming flux vector splitting, 25 which signif- 
icantly reduces the dissipation of the original scheme. 
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Time 

(**••) 

Altitude 

(km) 

R< D 

Velocity 

(krri/s) 

Density 

(kg/m 3 ) 

Too 

(K) 

1634 

76.42 

2x KV 1 

1 1 .36 

3. 72c- 5 

195 

i one 

71.04 

5 x I0' 1 

11.31 

8.57c-5 

210 

1637.5 

67.05 

7 x I0 4 

1 1 .25 

1.47c-4 

228 

163!) 

63.11 

1 x l() r ‘ 

11.14 

2.4 lc-4 

242 

16'IO.ii 

59.26 

2 x 10 r ’ 

10.97 

3.86c-4 

254 


Table 2 Fire II trajectory points and freest ream con- 
ditions. Rep is the freest ream Reynolds number based 
on body diameter, dime is given in seconds after 
launch. 

Third-order spatial accuracy is obtained through 
MUSCL extrapolation, coupled with aminmod limiter. 2e 
Time advancement to a steady state solution is achieved 
using the data-parallel line relaxation method. 24 This 
code has been previously applied successfully to several 
problems. 27 " 28 DPLR Axi typically required about 
10.000 iterations to reach a steady-state solution on 
the baseline (161 x 07) single block grid (as compared 
to about 1000 iterations for a forebody only case), 
and required approximately 33 CPU hours on a SGI 
R. 12000 workstation. The large scale flow features were 
set up fairly quickly, but the afterbody heat transfer, 
especially on the flat base plate, converged to its final 
answer very slowly. 

Wake flows can be sensitive to details of the volume 
grid construction, and care must be taken to generate 
a grid that is well aligned to the flow features. In par- 
ticular, it is important to accurately capture the rapid 
expansion around the shoulder, which determines the 
thermochemical state of the essentially frozen wake. 
The nature of the separated flow region is also de- 
termined in large part by the shear layer and wake 
compression region. 2 * At higher Reynolds numbers the 
wake will consist of mult iple counter-rotating vortices, 
which must be resolved. 2 * Finally, care must be taken 
to ensure that the grid completely encloses the sub- 
sonic portion of the wake, which can extend several 
body diameters downstream. For this paper, solutions 
were obtained using two grid topologies. The first 
topology, shown in Fig. 3 a, consists of a single mesh 
block constructed using a hyperbolic grid generator 31 
to ensure orthogonality at the surface. The solution 
adaptive grid code SAGE 32 is used to fit the bound 
ary to the computed shock. This topology allows for 
rapid grid generation and solution turnaround, and 
should be adequate for preliminary heat shield design, 
as well as for the computation of low Reynolds num- 
ber Hows where the wake structure is dominated by a 
single large vortex. However, as the Reynolds number 
increases and 1 he wake becomes more complex, it hs 


desirable to use a grid topology that allows for com- 
pact enrichment of the separated flow region without 
propagation of the grid points throughout the rest of 
the grid. The second grid (Fig. 3b), constructed us- 
ing GridPro 33 allows for compact enrichment with a 
topology that results in 36 elemental blocks, merge 
able into four. With such a topology the grid in the 
separation region can be arbitrarily densified without 
increasing the point density in the rest of the grid. 
With either topology, SAGE 32 is used to adapt the 
grid to the computed solution and cluster points near 
flow features of interest. 

Results 

Solutions were obtained for five trajectory points 
during the first data period, with freestream condi- 
tions listed in Table 2. The solutions range from 1634 
to 1640.5 seconds after launch, which spans the period 
from the first appreciable afterbody heat transfer mea- 
surements to a point near peak afterbody heating. So- 
lutions were limited to the first data period for several 
reasons, including the desire for laminar axisymmetric 
steady flow and the lack of complete geometric infor- 
mation for the second and third forebody calorimeter 
configurations. It should be noted that the final tra- 
jectory point, at t = 1640.5s, occurs after the onset of 
forebody calorimeter melting. Therefore, it is possible 
that the accuracy of the afterbody data may be com- 
promised by beryllium droplets entrained in the wake. 
However, the temperature vs. time traces for the after- 
body thermocouples show no anomalous readings that 
can be traced to forebody melting. 

Wake Structure 

Figure 4 shows streamlines in the separated flow re- 
gion for all five trajectory points. All solutions were 
obtained on the baseline single block grid, which has 
161 points along the body surface and 07 points in 
the normal direction. At the lowest Reynolds num- 
ber point (f — 1634s) the flow separates on the coni 
cal frustum well after the heatshield shoulder, and is 
dominated by a single large vortex. As the Reynolds 
number (and t) increases, the separated flow region 
becomes larger and more complex, as would be ex 
pected for a laminar flow. At f = 1636s and t — 1637.5s 
the flow separates at the rearward facing step after 
the shoulder, but remains dominated by a single vor- 
tex. However, by t- 1639s a bulge can be seen about 
a third of the way down the conical frustum, which 
can be further resolved into a small secondary vor- 
tex. Additional small vortices appear in the corners of 
the rearward facing steps at the rear shoulder and base 
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Figure 4. Streamlines in the separated flow region for 
solutions between t — 1 (3 3 '1 s and t — 1 64 0 showing 
evolution of the wake structure. 

of the probe. Finally, by t= 16d0.5.s a well developed 
secondary vortex is apparent on the frustum. While 
all computations predict these flows to be steady, it is 
likely that as the Reynolds number increases further 
the secondary vortices will begin to move and the flow 
will become unsteady.** 1 ' 

Surface Heating 

Figure 5a shows tiie computed total convective 
) and the catalytic component (q cat ) of heat trans 




b) 


Figure 5. Computed total convective (jtot) a nd cat- 
alytic component (g cat ) of heat transfer to the body 
surface at t — 1639s. a) Entire vehicle, and b) enlarge 
ment of the afterbody region with experimental data 
from the twelve calorimeters. 

fer to the surface at / — 1639s. Data are plotted vs. s, 
the pathlength distance from the apex of the vehicle. 
The surface was assumed to be non-catalytic to neutral 
species and fully catalytic to ions. The isothermal wall 
temperature (T w = 1325 K) was taken from the fore 
body thermocouple data. The solution was obtained 
on the baseline 161 x 97 grid, with the wall spacing 
chosen to maintain a constant cell Reynolds number 
(Rc c = 2) on the surface. The effects of grid resolution 
and varying wall spacing will be examined below. It is 
impossible to directly compare our forebody convec 
tive heating results to the experimental data, which 
are based on total (convective and radiative) heating. 
However, from Fig. 5a we see that the stagnation con 
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vective heating is about 330 W/cm 2 , which is within 
5% of the value reported by Olynick et.al. 4 using a 
similar set of modeling assumptions. The catalytic 
component of heating is negligible on the entire sur 
face, comprising about 2% of the total heat transfer. 
This is because most of the ions recombine in the cold 
boundary layer before reaching the surface. 

The flow undergoes a pair of rapid expansions at 
the heatshield shoulder (point 1 in Fig. 5a) and the 
rearward facing step (point 2), which result in a local 
increase in heat transfer during the expansion, followed 
by a rapid decrease. An enlargement of the afterbody 
region is shown in Fig. 5b. The numbers in Fig. 5b 
correspond to the inset in Fig. 5a. Also included in 
this figure are the experimental data from each of the 
12 calorimeters. Since ihe calorimeters were placed in 
a symmetrical array, there are three points at each s- 
location, corresponding to different circumferential (<j>) 
locations. The experimental data show that there is 
little circumferential data scatter, which is consistent 
with the small angle of attack (a < 1°) reported dur 
ing the first data period. However, the computation 
predicts a heat transfer rate that is considerably lower 
than the experimental data. Although not shown here, 
the same trend is observed for the other trajectory 
points. The rather poor comparison between the ini- 
tial computations and the data is not too surprising, 
considering the simplistic nature of the boundary con 
ditions that were initially imposed at the surface. In 
the next sections of the paper we look at several model 
ing assumptions that were made in order to determine 
their affect on computed heat transfer. 

Grid sensitivity 

One possibility for ihe poor agreement between the 
computation and experiment is grid resolution. In or 
der to examine this effect, grid refinement studies were 
conducted at the t~ 1634 and 1 = 1639 second trajec- 
tory points, covering a range of Reynolds numbers be- 
tween 2 x 10 4 to 1 x 10 n based on body diameter. The 
number of grid points in the axial and normal direc- 
tions as well as the body normal spacing were varied 
to determine their effect on computed heat transfer. 

Figure (> shows representative results at t = 1639s. 
Solutions are shown on three grids. Only the afterbody 
region is shown in Fig. 6; the forebody heat trans- 
fer varied by less than 1% on all grids tested (slightly 
larger deviations were seen at the stagnation point due 
to numerical issues with some of the grids). The base- 
line grid (“A’’) has lfil points along the body surface 
and 97 points in the normal direction. Wall spacing 
was chosen to maintain a constant cell Reynolds nmn 



Figure 6. Impact of grid resolution on the computed 
heat transfer at t — 1639s. The first grid is the baseline 
grid used in this paper. The next two solutions were 
obtained on single block grids with increasing axial 
resolution. The final case was obtained on a four block 
grid with compact enrichment of the wake core. 

ber ( Re c = 2), which implies that the near wall spacing 
in the low density afterbody is much larger than on the 
forebody. The second grid ( U B ) was refined in the ax 
ial direction only, with 253 axial points clustered near 
the shoulder and geometric corners on the afterbody. 
The normal spacing was identical to the previous grid. 
From Fig. 6 we see that the increased axial resolution 
lowered the predicted heat transfer by 10 — 15%. This 
trend is consistent with the results of Olynick,** 4 who 
showed that the corner expansion plays a significant 
role in determining heat transfer on the afterbody for 
the Stardust sample return capsule. Based on this re- 
sult, a third grid (“C”) was generated, which had 395 
points in the axial direction and the same normal dis- 
tribution as the previous grids. The solution on this 
grid was essentially the same as that obtained on grid 
B, indicating that grid B has sufficient axial resolu- 
tion. The final solution shown in Fig. 6 was obtained 
on a four block grid ( kt D v ) with 36,000 total grid points 
(50% more than grid B). Axial spacing was essentially 
identical to grid B (253 points on the surface); however 
the normal spacing was set to a constant value of 10 
m over the entire surface, and most of the additional 
points were clustered in the separated flow region using 
the compact enrichment technique described earlier. 
The computed heat transfer on this grid is essentially 
identical to that obtained on grid B, indicating that at 
these low Reynolds numbers the near wall spacing and 
normal resolution of the separated flow region are not 
critical to predicting heat transfer. No unsteadiness 
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Figure 7. Effect of the thermal state of the gas on 
computed afterbody heat transfer at i — 1639s. 

was observed in any of the solutions obtained for this 
study. 

A similar study was performed at t= 1634s. At this 
point, the Reynolds number is a factor of five lower 
than at t — 1639s\ and as expected there is less sensi- 
tivity to grid refinement. For this case the computed 
afterbody heat transfer varied less than 5% on all grids 
tested. These results indicate that grid resolution is 
not Ihe cause of the disagreement between the flight 
data and the computations. 

Sensitivity to the gas thermal state 

The effect of the assumed thermal state of the gas is 
shown in Fig. 7 for 1639s. The thermal nonequilib- 
rium solution is obtained using the method described 
above, while the equilibrium solution is computed us- 
ing statistical mechanics to model the vibrational en- 
ergy of the gas at the single temperature T. As seen in 
Fig. 7, the assumption of thermal equilibrium vs. non- 
equilibrium has no impact on the computed heat trans- 
fer. The same is true for all other trajectory points 
as well. This result is due to the fact that the post 
shook flow consists entirely of atomic species and their 
ions; there are essentially no molecules to carry the 
vibrational energy. Since the forebody wall boundary 
condition is non-catalylic to neutral species, there is 
little recombination occurring in the near- wall region 
on the forebody, and thus the wake core will consist 
mainly of atoms as well. For example, at / = 1639.*? 
the maximum mass fraction of diatomic species in the 
wake is less than 3%. However, for a less energetic 
flow, or for a forebody surface that is catalytic to 
neutrals, one would expect the thermal state of the 
gas to have a larger influence on the wake core. ,H 



Figure 8. Effect of the thermal wall boundary condi 
tion on computed afterbody heat transfer at / = 1639*. 

Sensitivity to surface temperature 

In the preliminary computations an isothermal wall 
was assumed, with the wall temperature taken to be 
the same as that for the forebody. However, the ma- 
terials used on the afterbody heat shield are much less 
conductive than the beryllium forebody. This would 
imply that perhaps a radiative equilibrium boundary 
condition would be more appropriate, in which con- 
duction into the surface is assumed to be zero and 
the heat converted to the surface is balanced by that 
radiated away. Heat transfer to the surface then be- 
comes a function of the surface temperature, given by 
q tot — c < tT? v . where e is the surface emissivity, rr i,s 
the Stefan-Boltzmann constant, and T ur is the wall 
temperature. At the other extreme, we can use the 
thermocouple data directly to predict an afterbody 
temperature. The thermocouple data indicate a tem- 
perature of about 400 K at /= 1639s. This tempera 
ture is slightly misleading, since the thermocouples are 
attached to highly conductive gold calorimeter slugs, 
rather than the much less conductive phenolic asbestos 
afterbody material. However, it provides another data 
point which can be used to assess the sensitivity of 
the heat transfer to the applied thermal boundary 
condition. Figure 8 shows the results of this analy- 
sis for t — 1 639s. Results are shown for a radiative 
equilibrium wall (e — 0.85), a “hot 5 * isothermal wall 
(T tu = 1325 K), and a cold wall (T w = 400 K). In each 
case the forebody temperature was held constant. The 
wall temperature on the conical frustum for the ra- 
diative equilibrium solution ranged from about 700 K 
at the beginning (s — 0.42 in) to about 1300 K at 
the end (s = 0.90 m). From the figure we see that 
decreasing wall temperature tends to increase surface 


i 





AIAA 200 1- 0065 9 



a) 



b) 

Figure 9. Assessment of continuum (low assump- 
tion, as measured by the gradient length local Knud- 
sen number (Kn^u). The continuum flow assump- 
tion is generally assumed to be fully valid in regions 
of the flow where Kw f n.i. < 0-05* a) £= 1634* and b) 

/ = 1639*. 

heat transfer, although the increase on the frustum is 
small (< 10%) as compared the the discrepancy be 
tween the flight data and the computation. Interest- 
ingly, the flat base (* > 0.91 m) was more sensitive to 
the wall temperature, with the heat transfer increas- 
ing by more than 20% as the wall temperature was 
lowered from its radiative equilibrium value (1500 K) 
to 400 K. However, since this region of the vehicle was 
simplified during construction of the OML, this effect 
was not explored further. 

Sensitivity to non-continuum effects 

The solutions shown in this paper were obtained as 
suming a continuum flow. However, at these Reynolds 
numbers non -continuum effects may be present in the 
base region. In order to assess the importance of non 
continuum effects on computed heat transfer, the gra- 


dient length local Knudsen number (An aiL ) 35 was 
computed for the t = 163d and t = 1639 second tra- 
jectory points, with the results shown in Fig. 9. Fol- 
lowing the work of Boyd et. al. 35 continuum break- 
down can be expected to begin when Kn aLL > 0.05. 
From Fig. 9a we see that at £ = 1634* [Reo = 2 x 10 4 ) 
there are significant non continuum effects in the near 
wall region of the afterbody. However, at £ = 1639s 
(Rep — 1 x 10 5 ) continuum breakdown is much less 
severe, and is confined to the last third of the con- 
ical frustum and the flat base (Fig. 9b). Unfortu- 
nately, Boyd et. al. did not quantify the relationship 
between I<n aLL and computed heat transfer. How- 
ever, a Navier-Stokes solution in general will slightly 
overpredict heat transfer in the base region of a non 
continuum flow. 36 Slip wall boundary conditions would 
permit a more rigorous estimate of non-continuum ef- 
fects with a Navier-Stokes code, and would capture 
this slight decrease in heat transfer. However, full im- 
plementation of a slip wall boundary condition in a 
reacting flow code is a non trivial effort, 3 ' which has 
not been completed in DPLR-Axi at this time. 

Sensitivity to surface catalysis 

The afterbody surface was initially assumed to be 
non-catalytic to neutral species. However, since the 
wake core consists almost entirely of dissociated atoms, 
any surface eatalycity would have a large impact on 
the heat transfer. Although the materials used in the 
construction of the Fire II afterbody heatshield are no 
longer in production, it is possible to draw some analo- 
gies to current materials. Based on analysis of silicon 
and carbon-phenolics, the surface coating and the un- 
derlying phenolic resin would likely start to pyrolyze 
when the surface temperature exceeded about 700 K. 
Vigorous pyrolysis would occur by about 800 K, and 
the process would be complete by the time the surface 
reached about 1100 K. 38 The end product of this py- 
rolysis would be a surface layer of carbon char, which 
should be nearly fully catalytic. However, this char 
would be mixed with asbestos fibers, which are essen 
tially a ceramic material with very low eatalycity. As 
stated previously, the surface temperature of the heat 
shield is not known accurately, making it is difficult to 
predict the onset of pyrolysis. However, the effects of 
surface eatalycity can be bounded by assuming a fully 
catalytic surface (7neutra/.t = 7 ion* = 1)* I R addition, a 
supposition 30 is made that, prior to the onset of pyrol- 
ysis, the surface eatalycity of the silicon based coating 
would be similar to that of SI RCA, 41 for which ac- 
commodation coefficients for nitrogen and oxygen re- 
combination as a function of temperature are known. 41 
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Figure 10. Comparison between heat transfer com- 
puted using several catalysis models on the afterbody 
with experimental data, a) / = 1634s and b) < = 1639s. 


Based on this analysis, solutions are obtained for 
each of the five trajectory points on the baseline grid, 
assuming an afterbody surface that is non-eatalytic, 
fully catalytic, and partially catalytic using the rates 
for SI RCA. In all cases a radiative equilibrium wall 
is assumed. Results for t = 1634s are shown, to- 
gether with the experimental data, in Fig. 10a. The 
assumption of catalysis on the afterbody significantly 
increases predicted heat transfer. As expected, 42 
both ( atalytic solutions show a pronounced catalytic 
“spike” at the rear of the shoulder (s — 0.38 rn), where 
the highly dissociated gas Hows over the surface tran- 
sition from a non catalytic to a catalytic surface. In 
addition, both solutions predict much higher heating 
on the flat base. This result is due to the fact that 
the flow is separated, and thus the highly dissociated 


wake flow first comes in contact with the catalytic sur 
face in the base. As the reverse flow travels forward 
on the frustum the flow rapidly recombines, deposit- 
ing its energy preferentially on the trailing end of the 
frustum. This effect is less pronounced for the SI RCA 
case, both because less catalysis is occurring and be- 
cause the SIRCA rates are temperature dependent. 
The hump in heat transfer for the fully catalytic so 
lution around s — 0.45 m is due to the fact the flow 
does not separate at the corner at this low Reynolds 
number (see Fig. 4), and thus there is a region of dis- 
sociated forward flow on the first part of the frustum. 
Once again, there is little scatter in the experimental 
data, with the exception of a spurious zero reading at 
$ = 0.63 m, <£ = 0° . The non-catalytic solution sig 
nificantly underpredicts the data, which is consistent 
with the results shown previously for tf=1639*\ How- 
ever, the fully catalytic solution slightly overpredicts 
the data, while the partially catalytic SIRCA solution 
is generally in good agreement. 

Figure 10b shows the same results for t = 1639s. For 
this case the catalytic jumps at the rear of the shoul- 
der and on the flat base are much larger, as expected. 
The SIRCA computation shows the best agreement 
with the data, although the heat transfer is overpre- 
dicted at the final two s-Iocations. The fully catalytic 
afterbody solution drastically overpredicts heating at 
these locations. Interestingly, the fully catalytic solu- 
tion predicts lower heat transfer on the forward part 
of the frustum for this case than the SIRCA solution. 
This effect appears to be due to the fact that for the 
fully catalytic case, much of the recombination (and 
associated heat release) is occurring at the back end 
of the frustum, while for SIRCA the slower recombi 
nation rates result in catalytic heat release along the 
entire surface. 

The trends in the data become more clear if we plot 
heat transfer vs. time at each calorimeter location, as 
shown in Fig. 11 for each of the four x/L calorimeter 
locations, where xj L — 0 at the rearward facing step 
on the heatshield shoulder and x/L — 1 at the theo 
retical apex of the conic frustum. In these plots, the 
solid lines are the experimental data from the three 
circumferentially distributed calorimeters at each x/L 
location, the dashed lines represent the assumed ±15% 
uncertainty in the experimental data, and the symbols 
are the results of the CFD computations for a non 
catalytic, fully catalytic, and SIRCA surface. From 
the figure, we see that the non-catalytic predictions 
are low throughout the trajectory. The fully catalytic 
and SIRCA predictions show good agreement with the 
data at x/L — 0.19 and xj L ~ 0.38 over the entire tra 
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Figure 11. Time histories of total heat transfer 
computed using several catalysis models on the after- 
body as compared to flight data, a) x/L = 0.19, b) 
x/L = 0.38, c) xj L = 0.56, and d) x/L- 0.70. x and L 
are defined in Fig. 1. 

jectory range examined here. However, at x/L — 0.58 
and x/L — 0.70 a different trend is seen. While the 
S1RCA prediction is in good agreement with the data 
up to t - 1637.5s, the fully catalytic solution has be- 
gun to significantly overpredict: heating by this time. 
In addition, as i increases, the heat transfer predicted 
by both catalytic solutions continues to rise, while the 
flight data appear to level off at 15-17 W / c.nr From 
a design standpoint, some comfort can be taken from 
the fact that the fully catalytic prediction is conser- 
vative over the later (highest heating) portion of the 
trajectory. However, the amount of conservatism in 
this solution on the last half of the frustum would lead 
to unnecessarily large margins on the heat shield thick- 
ness if it were used as a basis for design. 

Figure 12 shows the experimental heating data for 
x/L - 0.56 over the entire trajectory. Profiles at the 
other locations have a similar shape. Also shown in 
Fig. 12 is the total (convective and radiative) heat 
transfer measured at the forebody stagnation point. 
Gaps in the forebody heating data correspond to peri- 
ods of calorimeter melting and ejection. From Fig. 12 
it is clear that the afterbody heating on the frustum 
reached a maximum about six seconds earlier than 
the forebody stagnation point. This result is con- 
trary to computations used for Stardust heatslueld 
design, 34 and flight data from a thermocouple located 
on the conical frustum of the Mars Pathfinder entry 
vehicle, 43 which both showed that maximum heating 
on the afterbody frustum occurred concurrent with or 
after forebody peak heating. This discrepancy implies 
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Figure 12. Time histories of total heat transfer mea- 
sured during flight on the conical frustum at xjL — 
0.56, and at the forebody stagnation point. Gaps in 
the forebody heat transfer curve correspond to periods 
of heatshield melting and ejection. 

that there was some physical process occurring during 
during the Fire II entry that was not observed in the 
Stardust computations or the Pathfinder entry. It is 
difficult to identify the mechanism that caused the ex- 
perimental heat transfer to level off at around t = 1638s 
without a more detailed investigation of the material 
response of the surface. However, one possibility is 
that pyrolysis is beginning at around this time. The 
onset of pyrolytic blowing would reduce the heat trans- 
fer to the surface over that predicted here. Pyrolysis 
did not occur on the Pathfinder afterbody, since peak 
heat transfer rates were an order of magnitude smaller 
than those encountered by Fire II. Another possibil- 
ity is the assumption of a non-catalytic forebody. If 
the forebody were partially catalytic to neutrals, the 
concentration of atoms in the wake core would be re- 
timed , lowering the magnit ude of the catalytic jump at 
the trailing end of the frustum and on the base plate. 

Future Work 

This paper makes a start at analyzing the after- 
body data from the Fire II flight experiment with 
CFD. While good agreement was obtained over a por- 
tion of the trajectory assuming a catalytic afterbody, 
a great deal more work needs to be done. The prop- 
erties of the afterbody materials should be studied 
in detail, in order to improve the catalysis and ther- 
mal boundary conditions. Material response should 
be coupled to the CFD solutions to model the onset 
of pyrolysis. Three-dimensional simulations with cou- 
pled conduction would allow us to model the different 
materials that make up the calorimeter surface (oxi 


dized nickel-chrome) and the phenolic-asbestos after- 
body heat shield. Although the rearward facing ra- 
diometer did not show any incident radiation, there 
was a significant amount of radiative heating to the 
forebody. This radiation was shown to be weakly cou- 
pled to the flow, slightly altering the nature of the fore- 
body flowfield. 4 This coupling could affect the shoulder 
expansion process and thus the wake core. The rates 
used here for chemical and energy exchange processes 
were mainly derived from experimental data in corn 
preserve flows, and thus the universal applicability of 
these models in a base flow simulation, which combines 
a compressive and strongly expanding flow region, is 
uncertain. 2 " Finally, if data from later in the trajec- 
tory are to be examined, turbulence becomes an issue 
and appropriate turbulence models must be applied. 

Conclusions 

The wake flow and afterbody heating environment 
of the Project Fire 11 flight experiment were simulated 
using computational fluid dynamics. Results were ob 
tained over the first portion of the trajectory, spanning 
the period from the onset of significant afterbody heat- 
ing to the peak heating point. Computed heat transfer 
on the conical frustum portion of the afterbody was 
compared to the flight data. Although forebody heat- 
ing results showed excellent agreement with previous 
computations, the initial non-catalytic predictions of 
afterbody heating were about a factor of two below the 
experimental values. Grid resolution of the solutions 
was confirmed and some of the modeling assumptions 
used during the CFD simulations were examined as 
possible sources of error. This analysis pointed to the 
possibility that catalysis was occurring on the after- 
body heat shield. Computations were then made as- 
suming both a fully catalytic afterbody surface, and 
a partially catalytic surface using a material model 
for SI RCA, which was chosen to approximate the cat- 
alytic properties of the surface coating prior to the 
onset of pyrolysis. The computations including finite- 
rate catalysis, particularly the SI RCA model, are in 
good in good agreement with the data over the early 
portion of the trajectory. However, the catalytic re- 
sults are conservative near the peak afterbody heating 
point, especially on the rear portion of the conical frus- 
tum. Comparison of the computations and flight data 
from Fire II with previous data suggests that another 
mechanism, possibly pyrolysis, may be occurring dur 
ing this portion of the flight, reducing the total heat 
transfer over the predictions that do not include mate- 
rial response. Although more work needs to be done, 
these results give increased confidence in the ability 
of CFD to, both conservatively and accurately, pro 
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diet the afterbody heating environment of a planetary 
entry vehicle during the laminar portion of the trajec- 
tory. 
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